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ABSTRACT 

Eighteen RR Lyrae variables (RRLs) that he in the " 12M clump" identified 
by the QUEST survey have been observed spectroscopically to measure their 
radial velocities and metal abundances. Ten blue horizontal branch (BHB) stars 
identified by the Sloan Digital Sky Survey (SDSS) were added to this sample. 
Six of the 9 stars in the densest region of the clump have a mean radial velocity 
in the galactic rest frame (V^sr) of 99.8 and a = 17.3 km s~^, which is slightly 
smaller than the average error of the measurements. The whole sample contains 
8 RRLs and 5 BHB stars that have values of Vgsr suggesting membership in this 
stream. For 7 of these RRLs, the measurements of [Fe/H], which have an internal 
precision of 0.08 dex, yield ([Fe/H]) = —1.86 and cr = 0.40. These values suggest 
that the stream is a tidally disrupted dwarf spheroidal galaxy of low luminosity. 
Photometry from the database of the SDSS indicates that this stream covers at 
least 106 deg^ of the sky in the constellation Virgo. The name Virgo Stellar 
Stream is suggested. 

Subject headings: stars: variables: other — Galaxy: halo — Galaxy: kinematics 
— Galaxy: structure 
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1. INTRODUCTION 

A cornerstone of the popular hierarchical picture of galaxy formation is the growth of 
galaxies with time through multiple mergers. For example, simulations of this process in the 
framework of ACDM cosmology by Bullock & Johnston (2004), suggest that a large galaxy 
such as the Milky Way experienced ~100 mergers with dwarf satellite galaxies, and most 
of them occurred in the first few biUion years after the Big Bang. The advent of large- 
scale surveys of the galactic halo during the past decade has provided conclusive proof of 
co-moving groups of stars (stellar streams) exhibiting the kinematic and stellar properties 
of tidally destroyed dwarf galaxies. The most clear-cut of these are the streams from the 
Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy and the Monoceros stream (Newberg et al. 
2002; Ibata et al. 2003) but a few other less conspicuous, and more uncertain, have been 
reported. It is also clear that other galaxies, most notably M31 (Ibata et al. 2001), have 
experienced multiple mergers. The question is not whether mergers occurred, but whether 
or not the number of mergers and their properties match the predictions of the hierarchical 
picture. 

In this paper, we report spectroscopic observations of the "12M clump" of RRLs that 
was first revealed as an over-density in the QUEST RRL survey (Vivas et al. 2001; Vivas 
2002; Vivas & Zinn 2003; Zinn et al. 2004; Ivezic et al. 2005) and later in the Sloan Digital 
Sky Survey (SDSS) as an excess of F type main-sequence stars (Newberg et al. 2002). The 
12M clump is roughly centered at right ascension (RA) 12M (186°), spans the RA range of ~ 
175-200°, the -2.3° to 0.0° range in declination (Dec) of the QUEST survey, and the galactic 
longitude and latitude ranges of ~ 279-317° and ~ 60-63°, respectively. These RRLs have 
mean V magnitudes, corrected for extinction, equal to 16.9 ± 0.2, which with My = 0.55 
for RRLs, yields a distance from the Sun (r©) of ~19 kpc. Speculation on the origin of this 
clump has centered on the possibihty that it is part of the streams of stars from the Sgr 
dSph galaxy (Majewski et al. 2003; Martinez-Delgado et al. 2004) or the stream from a now 
extinct dSph galaxy (Vivas & Zinn 2003). 

2. OBSERVATIONS 

Our campaign to obtain spectroscopic observations of the RRLs in the 12^4 clump began 
in 2001, but only very modest success was obtained until 2005 because of poor weather and 
instrument malfunction. Eighteen RRLs have now been observed (QUEST # 177, 182, 189, 
191, 195, 196, 199, 200, 205, 209, 210, 224, 225, 226, 233, 236, 237, 261), and 12 have been 
observed more than once. With the ESO 1.5m and the Clay Magellan 6.5m telescopes, we 
obtained spectra centered at 4500A with resolutions of 3.lA and 4.3A, respectively. With 
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the WIYN telescope, we obtained spectra centered at 8550A with a resolution of 2.9A. 
Comparison lamp exposures were taken before and after each stellar one. Stars from Layden 
(1994) were also observed to serve as standard stars for radial velocity and for his pseudo- 
equivalent width system. With the blue spectra, we cross-correlated each target spectrum 
with ~10 spectra of several radial velocity standards of similar spectral type. With the red 
spectra, we fitted Gaussian line profiles to unblended Paschen lines of hydrogen. To measure 
the systemic velocity, we fitted the radial velocity curve of X Arietis to the type ab variables 
(see Vivas, Zinn, & Gallart 2005) and for the type c stars (191 & 224) we fitted a template 
that we constructed from velocity curves of T Sex and DH Peg. The errors in velocities 
were calculated taking into account the uncertainties in the phases of the observations, the 
template velocity curves, and the fits. These errors range from 8 to 22 km s~^, with a mean 
value of 17 km s^^. Because the RRLs are spread out over ~20° in RA, the Sun's peculiar 
motion and the motion of the Local Standard of Rest contribute different amounts to their 
radial velocities. We have therefore removed these effects and base the following analysis 
on the galactic standard of rest velocity (Vgsr), which is the radial velocity measured by an 
observer at the Sun who is at rest with the galactic center. To measure [Fe/H] from the blue 
spectra of the type ab variables, we followed closely Layden's variation of the Preston AS 
technique (Layden 1994), which cannot be apphed to red spectra or to type c stars. The 
average o"[kc/h] for stars observed more than once is 0.08 dex, which we consider the internal 
precision obtained with only one spectrogram (the external error is ~ ±0.15 dex, Layden 
1994). 

To this sample of RRLs, we added 10 BHB stars that were discovered by Sirko et al. 

(2004a) in the RA range 175-205°, the Dec range -2.5-0.0°, and the range 15-25 kpc. 
These stars occupy lines 22, 23, 28, 35, 38, 41, 112, 121, 126, and 774 in table 3 of Sirko 
et al. (2004a). We have adopted their values of Vgsr (Vios in their nomenclature) and also 
their values of Tq, which are on essentially the same distance scale as the RRLs. Sirko et al. 
(2004a) estimated that their velocities have an average a of about 26 km s^^. 



3. The Kinematic Signature of the Stream 

In Fig. 1, we have plotted our selection of RRLs and BHB stars together with the other 
QUEST RRLs in the region 175° <RA< 205° and 14< Tq <26 kpc, which encompasses the 
I2M clump. A very tight configuration of RRLs exists at RA ~186° and Tq ~ 19.6 kpc, 
and we have drawn a box around this "central region" in Fig. 1. We drew this box so as 
to include a reasonable sample (~ 1/3) of the total number of stars that had been observed 
spectroscopically and none that had not been observed. The model of the galactic halo that 
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Vivas & Zinn (2003) found was a good match to the distribution of type ab RRLs outside 
the densest regions in the QUEST catalogue predicts that this box should contain only 0.8 
of a type ab variable. Since the type c to ab ratio ranges from ~0.2 to ~0.8 in globular 
clusters, depending on Oosterhoff class (Smith 1995), the box is expected to contain ~1 to 
1.4 RRLs of all types. The observed number of 8 RRLs, 7 type ab and 1 type c, is clearly 
a large excess over the expected number. Expanding the box by either reducing the lower 
limit on by ~2 kpc or by increasing the upper RA limit by ~5° produces only a little 
dilution of the kinematic signature of the stream. 

In the upper panel of Fig. 2, we show the histogram of the V^gsr values of the 9 stars 
in the box in Fig. 1. The bins are 30 km s^^ wide, which is slightly larger than the la 
errors of the stars with least precise measurements. Also plotted is a normal distribution 
for (T^sr) = and a — 100 km s~^ . A random selection of halo stars is expected to a 
have a normal distribution (e.g., Harding et al. 2001) and values of (Igsr) and a close to 
these values (Sirko et al. 2004b; Brown et al. 2005). Harding ct al. (2001) recommended 
the Shapiro & Wilk (1965) statistical test of normality as a criterion to apply to halo fields 
that may harbor stellar streams. According to this test, the Vgsr values of the 9 stars are 
not normally distributed (>98% confidence), and this test is completely independent of any 
choices for the mean or a of the normal distribution. Six of the 9 stars have very closely 
clustered values of V^sr, {Vgsr) = 99-8, and remarkably a — 17.3 km s~^ , which is smaller 
than the average error (18.8 km s""*^) of the velocity for these stars. This tight velocity 
distribution in a volume of space where there is clearly an excess of RRLs is unequivocal 
evidence for a stellar stream. We suggest that it be called the Virgo Stellar Stream (VSS) 
after the constellation in which it is found. 

We do not expect to see such a narrow peak in V^gr as the sample is expanded because 
our measuring errors are substantial, the stream probably has some intrinsic dispersion in 
velocity, and perhaps also a gradient in (V^sr) with position. Therefore to investigate the size 
of the stream, we consider as possible member any star whose V^sr lies within ±60 km s~^ 
of the (V^gsr) of the central region (i.e., 40 < Vgsr < 160 km s^^, between the vertical dashed 
lines in Fig. 2). In terms of the measuring errors, this range corresponds to ±2.3(T and ±3.6cr 
for the BHB stars and RRLs, respectively. Consequently, in absence of a sizable velocity 
dispersion and/or gradient, this range is likely to include all the members. 

The histogram of the values of Vgsr for the stars lying outside of the central region is 
shown in the lower panel of Fig. 2. The kinematic signature is much diluted, but it is not 
completely absent, for there is an excess of stars in V^sr range expected of members. Several 
of these stars have lower values of Vgsr than the members in the central region. This could 
be due to a velocity gradient, but it is also consistent with the measuring errors and the 
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contamination by non-members,which are more likely to have small values of T^sr than larger 
ones. We have plotted as solid symbols in Fig. 1 the 13 stars^ in the whole sample that meet 
the Vgsi criterion for membership. These stars span large ranges in RA and Tq, suggesting 
that the VSS is both large and diffuse. Two of these stars stand out as the most likely 
non-members: BHB star 23, which lies ~4 kpc more distant than the major concentration, 
and RRL 261, which is the RRL lying farthest to the East and is also the star in this sample 
with the largest value of Vgsr ■ The removal of these two stars from the sample changes (V^sr) 
and a from 85.4 and 31.4 km s~^, respectively, to 83.2 and 24.5 km s~^. Since this last a 
is not much larger than the average measuring error for this sample (20.3 km s~^), more 
precise measurements are needed to refine the membership criterion and to be certain that 
the velocity dispersion has been resolved. Our future measurements of RRLs should indicate 
whether the VSS extends in the directions suggested by BHB 23 and RRL 261. 



4. The Signature in the Luminosity Function 

Newberg et al. (2002) discovered an excess of F type main-sequence stars in SDSS in 
the direction to the 12^4 clump and at the appropriate magnitude to be related to the 
RRLs (feature S297+63-20.0). In fig. 1 of Newberg et al. (2002), the globular cluster Pal 5 
appears as a streak pointing radially toward the Sun. This is a sign that this method, which is 
excellent for finding halo substructure, may not locate precisely an old stellar population. We 
have therefore examined the SDSS photometry for evidence of the upturn in the luminosity 
function (LF) of an old stellar population that occurs near the magnitude where the subgiant 
branch (SGB) and the main-sequence (MS) merge (see also Martinez-Delgado et al. 2002). 

The method is illustrated in the top 3 panels of Fig. 3. The top left panel shows the 
color magnitude diagram (CMD) of 2°x2° centered on the globular cluster Pal 5. The cluster 
and its tidal tails occupy only a fraction of this field. The top middle panel is the CMD of 
a control field of equal area. The top right panel shows the LFs of the Pal 5 field and the 
control field. These functions were constructed using the color indices cl and c2 as defined by 
Odenkirchen et al. (2001) and their filtering technique, which minimizes the contamination 
by poor measurements and by field stars. While the presence of Pal 5 is clearly detectable 
in the CMD itself, it is very evident in the right panel as a sharp increase in the number 
of stars at g~20.25, which is roughly the magnitude of the SGB (the MS turn-off is ~ 0.5 
fainter. Smith et al. 1986). 



^Eight of them are RRL stars (177, 189, 195, 196, 199, 210, 237 & 261) and 5 are BHB stars (# 22, 23, 
28, 38, 121, in table 3 Sirko et al. 2004a) 
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The middle panels in Fig. 3 show the CMDs and the LFs of the central part of the 
VSS and a comparison field. The 3 RRLs in this VSS field have {V) = 17.09, whereas the 
5 RRLs in Pal 5 have {V) = 17.44 (Vivas et al. 2004). If VSS resembles Pal 5, a sharp 
upturn in the LF is expected at g~19.9. While the VSS field shows an excess of stars over 
its comparison field at this magnitude, its LF clearly does not mimic that of Pal 5. However, 
Pal 5 is a moderately metal- rich globular cluster ([Fe/H] = -1.47, Zinn & West 1984), and 
it has the combination of age and metallicity that produces a nearly horizontal SGB and 
consequently a very sharp upturn in the LF . The 7 RRLs in the VSS with measurements of 
[Fe/H] have a mean value of -1.86, which suggests that a more metal-poor globular cluster 
than Pal 5 would be a better comparison object. The LF constructed by Zaggia, Piotto, & 
Capaccioh (1997) for the globular cluster M55 ([Fe/H] = -1.82, Zinn & West 1984) shows a 
more gradual increase because its SGB is more steeply inclined in the CMD than Pal 5's. 
If the LF of VSS resembles that of M55, a relatively gradual climb in the number of stars 
is expected to start at My ~ +3.0, or g~ 19.5. Since the LF in the middle panel of Fig. 3 
shows this behavior, we believe the excess of F type main-sequence stars first recognized in 
the SDSS data by Newberg et al. (2002) is indeed the MS of the VSS. The bottom panels 
show the CMDs and the LFs for a field 10° to the East of the central one, where the VSS 
appears to be present but not as strongly as in the central one. 

Using this technique with the SDSS photometry, we searched for the VSS in 272 deg^, 
between 176° < RA < 210° & -4° < Dec < +4°. We considered that the VSS was detected if 
in the range 19.5 < g < 20.5 the number of stars in the target field was consistently greater 
than the number in the control field by an amount larger than the combined Poisson errors. 
The control and target fields had similar galactic latitudes but were offset from each other 
by 10 or more degrees. Fig. 4 shows the region of the sky (~ 106 deg^) where the VSS was 
detected. Towards the east, the VSS becomes less prominent due to the progressive increase 
of the field population. The northen part of the VSS (0° < Dec < 2°) has a somewhat 
brighter upturn in the LF than the southern portion (—4° < Dec < —2°), which indicates 
that it is closer to the Sun. 



5. Discussion 

The previous suggestions (Majewski et al. 2003; Martinez-Delgado et al. 2004) that the 
I2M clump of RRLs may be part of the streams from the Sgr dSph are inconsistent with our 
observations. The model of the Sgr streams by Martinez-Delgado et al. (2004) predicts a very 
low density of Sgr stars in this volume of space, which is in conflict with the concentration 
seen in Fig. 1. The models by Law, Johnston, & Majewski (2005) that assume either a 
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spherical or an oblate shape for the Milky Way's dark matter halo, predict more significant 
numbers of Sgr stars in this volume but with (V^sr) ~ -180 and ~ -260 km , respectively. 
These values are completely incompatible with the (V^sr) of the VSS (+83 or +100 km 
, depending on sample selection). According to the prolate model of Law, Johnston, & 
Majewski (2005), this region should not contain debris from Sgr. 

Based on our measurements for 7 type ab RRLs that have values of V^sr that are consis- 
tent with membership, the VSS has ([Fe/H]) = —1.86 and a = 0.40. Because this dispersion 
is several times the average of the o"[Fc/h] values (0.08), we conclude that the VSS has a sig- 
nificant range in [Fe/H]. A range of this magnitude is characteristic of all the dSph satellite 
galaxies of the Milky Way, but not of the vast majority of globular clusters. The VSS is 
probably the debris of a disrupted dSph galaxy. There is no sign of a significant interme- 
diate age population in the CMDs from the SDSS data (Newberg et al. 2002 and the ones 
we produced), and the fists of halo carbon stars (Totten & Irwin 1998; Mauron et al. 2004; 
Mauron, Kendall, & Gigoyan 2005) do not contain any that are consistent with membership 
in the VSS. Thus, the progenitor of the VSS appears to have been a system dominated by its 
very old stellar population. These properties suggest that of the extant dSph galaxies, the 
VSS may most closely resemble the ones of the lowest luminosity, e.g., Ursa Minor, Draco, 
and Sextans (Mateo 1998). 

Finally, our measurements have also revealed a few other interesting features. Ignoring 
the VSS stars, there arc some smaller groups of stars that have positions and values of (V^gsr) 
that suggest possible membership in the same "moving group" (e.g., RRL 226, 233, & BHB 
35; RRL 224, 225, 236 & BHB 126). Furthermore, if we consider together the 15 stars that 
are not members of the VSS and the two that have low probabifities of membership (RRL 
261 & BHB 23), we find (^sr) ~ —9.1 and a — 163.7 km s~^ . This value of a is remarkably 
large, and it is inconsistent at above the 95% confidence level (F test) with the value of 101.6 
km that Sirko et al. (2004b) found from their sample of 1170 BHB stars with median 
~ 25 kpc. The origin of this kinematically hot distribution, the possibility of other streams 
besides the VSS, and the relationship of the VSS to the large feature recently reported by 
Juric et al. (2005) in Virgo are being investigated as we obtain additional spectroscopy of 
QUEST RRLs. 

This research is part of a joint project between the Universidad dc Chile and Yale 
University, which is partially funded by the Fundacion Andes. We acknowledge funding from 
the following institutions and grants: NSF AST 00-98428, NSF AST 05-07364, Fundacion 
Andes C-13798, MECESUP UCH-118 and FONDAP No. 15010003. 
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Fig. 1. — The distance from the Sun (r©) is plotted against right ascension (RA) for the 
RRLs in the QUEST survey (circles) and for the BHB stars (triangles) discovered by Sirko 
et al. (2004a). Large circles depict the RRLs that were observed spectroscopically. Sohd 
symbols depict the RRLs and BHB stars that are probable members of the Virgo Stellar 
Stream. The dashed lines enclose the central region (see text). The errors in scale with 
distance, and the plotted error bar is appropriate for RRLs at 20 kpc. The errors for the 
BHB stars are roughly a factor of 1.4 larger. 
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Fig. 2. — The histograms of the Vgsr vahies of the stars inside and outside the central region 
are plotted in the top and bottom panels, respectively. Random samples of halo stars that 
contain the same number of stars are expected to have the normal distributions (solid lines). 
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Fig. 3. — In the left panels, we show the CMDs from SDSS photometry for 2° x 2° fields 
that include (top) the globular cluster Palomar 5, (middle) the center of the VSS (186° < 
RA < 188°, -2° < Dec< 0°), and (bottom) an eastern field of the VSS (196° < RA < 198°, 
-2° < Dec < 0°). In the middle panels, we show the CMDs for control fields of equal size 
that were offset by 10 or more degrees. In the right panels, the luminosity functions of the 
object field (solid line) and the control field (dashed line) are plotted as histograms (Poisson 
error bars). 
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Fig. 4. — The dashed hnes enclose the region of the sky where we found an excess of stars 
compared to the control fields in the interval ~ 19.5 < g < 20.5. The RRLs and BHB stars 
that are possibly members of the Virgo Stellar Stream (i.e. the stars within the dashed 
vertical lines in Fig. 2, and solid symbols in Fig. 1) are plotted as circles and triangles, 
respectively. 



